Electron-phonon Interaction close to a Mott transition 
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The effect of Holstein electron-phonon interaction on a Hubbard model close to a Mott-Hubbard 
transition at half-filling is investigated by means of Dynamical Mean-Field Theory. We observe a 
reduction of the effective mass that we interpret in terms of a reduced effective repulsion. When 
the repulsion is rescaled to take into account this effect, the quasiparticle low-energy features are 
unaffected by the electron-phonon interaction. Phonon features are only observed within the high- 
energy Hubbard bands. The lack of electron-phonon fingerprints in the quasiparticle physics can be 
explained interpreting the quasiparticle motion in terms of rare fast processes. 
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Electron-phonon (e-ph) interaction is a key effect 
in solids, that gives rise to superconductivity, charge- 
density-waves and other important effects. While the 
physical phenomena caused by e-ph coupling in weakly 
correlated materials are reasonably well understood, 
much less is known about the role of e-ph coupling in 
strongly correlated materials, i.e., in compounds in which 
the short-range Coulomb repulsion is larger than the typ- 
ical kinetic energy. This issue, besides its obvious gen- 
eral interest, has been strongly revived by recent experi- 
ments which suggest a role of lattice effects in the high- 
temperature superconducting cuprates (which are known 
to be strongly correlated compounds) I'l , and by the real- 
ization that strong correlations are important inphonon- 
mediated superconductors as the fullerenes 0, 0] . 

The interplay between e-ph and electron-electron (e-e) 
interactions is quite a complicated problem, in which a 
variety of physical regimes and different phases can be 
realized. In the present paper we focus on the strongly 
correlated regime, in which the e-e interaction is assumed 
to be quite large, and the mobility of the electrons is 
strongly reduced. Our main concern is to understand 
whether and to which extent the e-ph coupling can still 
influence the electronic properties in this regime. We no- 
tice that the symmetry of the e-ph coupling is important 
in this regard: It has indeed been shown that a Jahn- 
Teller coupling, where the phonon variables are coupled 
with orbital and spin degrees of freedom is weakly af- 
fected by Hubbard repulsion [/, whereas a Holstein cou- 
pling with the local charge is strongly suppressed 0, . 
Nevertheless, even the Holstein coupling, if sizeable, can 
have a qualitative physical effect in a correlated material. 
For instance, it can induce phase separation close to the 
density-driven Mott transition A recent important 
contribution comes from the two-dimensional Quantum 
Monte Carlo study of Ref . 0] , where a small U is found to 
suppresses the e-ph coupling at all electron and phonon 
momenta, while, in the strongly correlated regime, the 
forward-scattering amplitude due to e-ph processes sub- 



stantially increases with U, leading to a non-monotonic 
behavior of the effective coupling. 

In this work we consider again the simplest case of the 
Holstein interaction, whose formal simplicity allows for 
extensive systematic studies. Our Hamiltonian reads 

-gY.^Ma■^+al)+UJoY,^ala„ (1) 

where Ci^cr (cj^.) ^'^d (a|) are, respectively, destruc- 
tion (creation) operators for fermions with spin a and 
for local vibrations of frequency t^o on site i, t is the hop- 
ping amplitude, U is the local Hubbard repulsion and g 
is an e-ph coupling. Despite its simplified character, the 
Hubbard-Holstein model is expected to display a com- 
plicated phase diagram as a function of its parameters, 
the correlation strength U/D {D is the electronic half- 
bandwidth), the e-ph coupling A = JoVwoi*, the adia- 
batic ratio loq/D and the density n [5|,l7j,|3- As notable 
examples, this model presents the physics of the Mott 
transition when the Hubbard repulsion is the most impor- 
tant scale ^ , a crossover to polaronic carriers or a bipo- 
laronic metal-insulator transition for small U and large 
A I?!, 01 , and it is expected to be unstable toward su- 
perconductivity, charge-density waves and antiferromag- 
netism in different regions of the phase diagrams 11]. To 
gain physical insight into some aspects of this phase di- 
agram, it is necessary to fix some parameters, and put 
ourselves in a given regime. As we have said, we investi- 
gate the strongly correlated metallic region, close to the 
Mott transition at half-filling. Our analysis is therefore 
relevant for metallic phases of correlated materials, such 
as the doped cuprates. 

We solve the model by means of Dynamical Mean-Field 
Theory (DMFT), which has emerged as one of the most 
reliable tools from the analysis of both correlated materi- 
als, and e-ph interactions. The method maps the lattice 
model onto a local theory which still retains full quantum 
dynamics, and it is therefore expected to be quite accu- 
rate for models with local interactions, such as 1^. The 
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DMFT solution is enforced by solving the local dynamical 
theory and imposing a self-consistency condition 9] . The 
local problem is equivalent to an Anderson-Holstein im- 
purity modelH^I which has to be solved with some tool. 
In our work we use exact diagonalization i truncating 
the infinite phonon Hilbert space allowing up to N„iax 
phonon states (ranging from 20 to 40), and using A^f, = 9 
sites in the conduction bath. 

DMFT has been extensively used to study the half- 
filled Hubbard model, showing that, in the metallic re- 
gion close to the Mott transition, the spectral function 
presents three features: the high-energy Hubbard bands 
at energy ±/7/2 and a narrow resonance at the Fermi 
level, whose width is proportional to the quasiparticle 
residue Z ~ (1 — dl2{u}) / du!\uj=o)~^ , a quantity which 
decreases for increasing U and vanishes at the Mott 
transition. Due to the momentum independence of the 
self-energy, Z is also the inverse of the effective mass 
Z = m/m*. On general grounds, we can expect two 
main effects of e-ph coupling starting from the situation 
described above. The first is a reduction of Z (equiv- 
alent to an enhancement of the effective mass). This 
effect is responsible for the polaron crossover for weak 
repulsion. The second effect is a phonon-mediated re- 
tarded attraction between electrons with opposite spins, 
which opposes the Hubbard repulsion leading to an ef- 
fective dynamical interaction Ueff{uj) = C/ — J^"2 . In 
the antiadiabatic limit wq — > oo the attraction becomes 
static and equals the bipolaronic binding energy 2g'^ /luq. 
While the first effect opposes to the electronic motion, 
the second effect instead favors the motion by screening 
the Hubbard repulsion and leads to an enhancement of 
Z (at least as far as C/ > 2g'^/ujo). The overall effect 
of the e-ph coupling on the correlated metallic phase is 
therefore hard to predict on intuitive grounds. 
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the quasiparticle weight Z as a function of A for different 
values oiU/D for the half-filled Hubbard-Holstein model 
in the relatively adiabatic regime ujq/D = 0.2. In the 
absence of Hubbard repulsion, the effect of the e-ph cou- 
pling is quite naturally to increase the electron effective 
mass, and decrease Z, eventually reaching the bipola- 
ronic metal-insulator transition [lOj . Then we turn on the 
Hubbard repulsion, and plot the ratio Z{U, X)/Z{U,0), 
in order to disentangle the effect of the e-ph interaction. 
While Z{U,X) is a decreasing function of U for each A, 
the above ratio displays a richer behavior. Increasing 
the value of U, the effect of the e-ph interaction becomes 
weaker, until a value of J7 ~ 1.5D is reached, for which 
the e-ph coupling has the surprising effect to increase Z 
(reduce the effective mass), in agreement with the re- 
sults found with a Numerical Renormalization Group so- 
lution of DMFT in RefiTj. This naively surprising re- 
sult can be attributed, in light of the previous discus- 
sion, to the reduction of the effective repulsion, which 
prevails over the polaronic renormalization of the hop- 
ping. In the following we empirically determine an ef- 
fective static repulsion Uef f(U, XjLOq), as the value of U 
for which a pure Hubbard model has the same Z as the 
full Hubbard-Holstein model. Since in the antiadiabatic 
limit Uef f = U — XD, we try to parameterize the effective 
repulsion as Ueff = U — rjXD, where 77 is a dimension- 
less free parameter that we determine by fitting directly 
Z{U, A) with this functional form. 
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FIG. 2: (Color online) The coefRcient 77 which measures the 
effective reduction of the Hubbard repulsion according to 
Ueff — U — rjXD, plotted as a function of uio/U. Different 
symbols refer to different values of U. The values of A used 
for the fitting procedure described in the text, range typically 
from to 1. 



FIG. 1: (Color online) Effect of e-ph interaction on the quasi- 
particle weight Z. The ratio Z{X,U)/Z{0,U) is shown for 
different values of U and uq = 0.2D. 

The DMFT study of model Q helps us in understand- 
ing which of the two effects prevails. In Fig. ^we report 



Computing 77 for different values of U and loq, we find 
the behavior shown in Fig. |21 For relatively small ojq/U, 
we find that 77 has a linear behavior in ujq/U, and bends 
down for larger phonon frequency closely following the 
functional form 77 = 2a;o/C//(l -I- 2ujo/U) and eventually 
reaching the asymptotic value 77 = 1. We will see in the 



following that this functional form can be derived from 
an expansion of the effective Kondo coupling [l^. 

The analysis of the quasiparticle weight suggests that 
the effect of e-ph interaction on the strongly correlated 
metallic phase of the Hubbard model is a partial screen- 
ing of the repulsion, and the degree of screening is con- 
trolled by uJa/U when this quantity is not too large. This 
latter factor contains the information about the lattice 
dynamics in an effective way. Of course, at this level, this 
result could sound of little physical significance, since it 
basically arises by a simple fit procedure. Nonetheless, Z 
also measures the width and weight of the low-energy 
feature in the single-particle spectral function. Thus, 
it is tempting to compare the spectral function for the 
Holstein model with the correspondent quantity for the 
equivalent Hubbard model with reduced repulsion. As 
far as [/ ^ tJoj we find that the low-energy part of the 
spectrum is basically identical in the two models, while 
differences develop in the high-energy part, where phonon 
satellites at energies of the order appear in the system 
coupled to phonons. This can be seen in FigI3 where the 
spectra for three different cases (see caption) are shown. 
Moreover, even if the effect of phonons is visible in the 
high-energy Hubbard bands, their position is more or less 
coincident in the Hubbard-Holstein model and in the ef- 
fective Hubbard model. This is quite surprising, since 
the rescaling of U has been derived only be requiring the 
same value for the quasiparticle weight. 

This findings strongly suggest that our reduction of the 
dynamical effective interaction to a static quantity is not 
a mere fitting, but it unveils some basic physics. We can 
gain some understanding of this scenario by exploiting 
the fact that the quasiparticle peak is associated within 
DMFT to a Kondo effect in the impurity model. It has 
been recently pointed out that, at the leading order in 
I/C/, the Kondo coupling Jk is only weakly affected by 
e-ph coupling 1^ . This result is closely related to a pre- 
vious calculation addressing the effect of e-ph interaction 
on the superexchange coupling in a Hubbard model . 
In both cases, the coupling (Kondo or antiferromagnetic) 
arises in fact, at the leading order in 1/U , from virtual 
processes in which doubly occupied sites are created, with 
and energy cost of order U . In the presence of an e-ph 
coupling, the intermediate states can have an arbitrary 
number of phonons, but, in the limit in which «C U , 
the corrections are small, and Jk is not substantially 
affected. More specifically, summing the contributions 
from all the states with different phonon occupation (see, 
e.g., Eq.(8) of RefEi), one obtains a small correction of 
Jk given, at first order in Xq/U, by 
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FIG. 3: (Color online) The electron spectral function for the 
Hubbard-Holstein model (dashed line) , compared with a pure 
Hubbard model with U — Ueff — U — rjXD (solid line), for 
different values of U and u>o (the dot-dashed vertical line in- 
dicates luq). In the two upper panels we take A = 1, in the 
lowest A — 0.5. 



which at first order in luq/U becomes Jk{X)/ Jk{0) — 
1 + 2ujqXD/U'^ ^ I + 4gVC^^- Eq. © determines 



an effective static U for the Kondo problem given by 
U^j:f = U~ifXD, with 77^ = {2uoo/U)/{l+2uJo/U). The 
agreement between the 77 extracted from DMFT and this 
prediction is evident from Fig. [3 This provides a strong 
physical interpretation of our previous finding about the 
effective U giving rise to the renormalization of the quasi- 
particle weight for small Lo^jV , for which we obtain the 
same dependence (see Fig. Indeed within DMFT, 
the formation of the quasiparticle peak can be associ- 
ated to a Kondo effect, which may in turn be described 
as arising from virtual processes of the same kind as the 
ones that lead to the Kondo coupling. But for the overall 
rescaling of Vef f , looking in detail to the phonon-induced 
modifications of the low-energy spectrum (the Kondo res- 
onance in the single-impurity language) we found that 
these modifications are quantitatively extremely small 
and almost invisible on the resolution scale of Fig|31 I(|. 

The picture that emerges from the previous arguments 
can be summarized as follows: Quasiparticle motion 
arises from virtual processes in which doubly occupied 
sites are created. Obviously, these processes are not so 
frequent, since the energy scale involved is large, but they 
are extremely rapid (the associated time scale is cx l/U), 
and consequently are poorly affected by phonon excita- 
tions with a characteristic time scale I/wq 3> l/U. When 
the phonon frequency is small with respect to U, the 
phonon degrees are frozen during the virtual excitation 
processes. Therefore, despite the overall electron motion 
is quite slow due to the small number of virtual processes 
(which is reflected by the large effective mass), the e-ph 
interaction has no major effect except for a slight reduc- 
tion of the total static repulsion. 

Finally, we come to some recent experimental evi- 
dences of phononic effects in the strongly correlated 
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cuprate superconductors. It must be noted that our re- 
sults refer to the half-filled case for U smaller than the 
Mott transition point, while the experiments obviously 
refer to a doped system most likely with a larger U. 
Nevertheless it is tempting to assume that also in the 
U vs. doping region, relevant for the cuprates, the main 
low-energy effect of phonons is a small down-shift of the 
repulsion, that will be of no relevance when U is appre- 
ciably larger than its value at the Mott transition point. 
Our main finding of a completely different effect of e-ph 
coupling on the low-energy part of the spectrum with 
respect to the high-energy region, may provide an ex- 
planation of the observation of e-ph interaction in the 
high-energy branch of the ARPES spectra, in contrast 
to the weak effect on the low-energy part, as observed 
in recent isotope effect measurements [l7|. On the other 
hand, if the system is close to the Mott transition, a slight 
modification of Ue/f = U — 2ujo\D/U through the vari- 
ation of loq can significantly change the effective mass, 
leading to important isotope effects, affecting the pene- 
tration depth, as found in various underdoped samples 
[lil . A more detailed analysis of the evolution of the iso- 
tope effects in the full U vs. doping diagram is worth 
being carried out to provide a more reliable comparison 
with the available experiments. 

A main limitation of our work, in light of previous 
studies of the problem of e-ph interactions, is the in- 
ability of DMFT to capture the momentum dependence 
of the phonon vertex. Even for the Holstein model, in 
which the bare coupling is local, i.e., momentum inde- 
pendent, it has indeed been shown that correlation ef- 
fects tend to favor small exchanged momenta (forward) 
scattering In this perspective, our DMFT results 

indicating that the local effects of the Holstein coupling 
are strongly renormalized down by correlations, allow to 
conclude that "standard" features of the e-ph coupling 
(polaron crossover, ordinary BCS pairing) are negatively 
influenced by correlation, while non-standard effects re- 
lated to the strong forward scattering can survive and 
even get enhanced. In other words, if any Holstein- 
like e-ph interaction turns out to be relevant in strongly 
correlated materials, it must display anomalous features 
associated to the relevance of forward scattering and 
consequently requires a wider framework than ordinary 
Migdal-Eliashberg theory. On the other hand, the strong 
reduction by U does not apply to e-ph coupling involv- 
ing degrees of freedom non competing with local charge 
fluctuations. This is the case for the Jahn- Teller coupling 
considered in Rcf. "s"] , where the phonon- mediated pairing 
is indeed enhanced by strong correlations. 

We have studied the half-filled Hubbard-Holstein 
model in the strongly correlated regime, where correla- 
tion effects are most important and determine a clear 
separation between low-energy quasiparticle features and 
"insulating" Hubbard bands. In this regime, we find 
that e-ph interaction has basically no effect on the quasi- 



particle features as long as the Hubbard repulsion is 
suitably rescaled. An effective static repulsion Ueff = 
U — 2loq\D/U in fact reproduces the low-energy features 
of the Hubbard-Holstein model using a simple Hubbard 
model. On the other hand, the high-energy Hubbard 
bands display phononic features that can not be found 
in an effective purely electronic model, even if the same 
scaling of U also reproduces the correct position of the 
bands. The "protection" of the low-energy physics with 
respect to e-ph interactions is provided by the (small) ef- 
fect of lattice coupling on the Kondo coupling and on the 
resonance shape. We understand this result in terms of 
freezing of phonons on the time-scale relevant to virtual 
processes underlying the quasiparticle motion. 
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